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In  a  study  of  perforated  muzzle  brakes,  Nagamatsu,  Choi,  Duffy,  and  Carofano 
calculated  the  three-dimensional  steady  flow  through  one  vent  hole  and  used 
the  results  to  predict  overall  brake  performance.  In  the  present  study,  the 
analysis  is  extended  to  the  calculation  of  the  blast  field.  The  results 
compare  favoraU3ly  with  previously  unpublished  shadowgraphs  obtained  by 
Dillon  in  his  experimental  program.  i\< vi.. fv'i' 
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Figure  1.  Schawatlc  drawing  of  a  perforated  nuzzle  brake. 

Nagaaatsu,  Duffy,  Choi,  and  Carofano  (raf  6)  calculated  tha  steady  three- 
dlaensional  flow  through  a  single  vent  In  tha  wall  of  a  shock  tuba.  Tha  pre¬ 
dicted  pressure  distribution  on  tha  vent  wall  conparad  favorably  with  tha 
axperiaental  aaasuraaants  of  Nagaaatsu,  Duffy,  and  Choi  (raf  7).  It  was  also 
shown  that  these  results  could  ba  eoabinad  with  a  ona-diaansional  aodal  of  the 
transient  flow  in  a  cannon  to  predict  tha  iapulse  reduction  produced  by  a  per¬ 
forated  muzzle  brake.  The  predictions  agreed  we11  with  the  expariaantal 
aaasi' *eaents  of  Di11on  (ref  1)  for  a  20-na  cannon.  A  more  coaprahansive  coa- 
parison  of  the  theory  with  these  expariaants,  including  a  discussion  of  brake 
efficiency  and  scaling,  was  aade  by  Carofano  (ref  8). 

References  are  listed  at  tha  end  of  this  report. 
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Th«  tr«ns1«nt  dtviopwtnt  of  tho  throo-dlMns^onol  flOM  thr.High  •  eonvontlonal 
bofflo  brokt  Is  diseustod  by  Msng,  Midhof>f«  and  Chon  (rof  9). 

In  tho  prosont  study,  tho  analysis  is  oxtondod  to  tho  calculation  of  tho 
blast  fio1d.  Tho  cooplox  transiont  throo*di»onsiona1  probloo  is  roducod  to  a 
ono*disonsiona1  floH  insido  tho  tubo  couplad  to  an  axisysnotric  oodol  outsido. 
Tho  throo-dioonsional  character  of  tno  flON  through  tho  tubo  wall  is  rotainod, 
howovor,  by  treating  it  as  quasi-stoady.  Tho  results  are  cooparod  with  pre¬ 
viously  unpublished  shadowgraphs  obtbinod  by  Dillon  in  his  oxporisontal  progr^o. 

THE  VENT  FLON  FIELD 

When  tho  propellant  gas  expands  through  tho  brake,  an  asyoootric  pressure 
distribution  develops  in  each  hole  with  the  highest  pressures  acting  on  the 
downstreas  surface.  The  vector  and  pressure  contour  plots  of  Figures  2  and  3 
show  typical  flow  patterns  in  the  synsetry  plane  of  one  hole  and  the  portion  of 
tho  tubo  associated  with  it.  Tho  flow  variables  in  tho  tube  are  uniforn  across 
tho  entrance  plane.  Tho  solid  linos  in  tho  vector  plot  indicate  where  the  local 
Mach  nuobor  is  unity. 

In  Figure  2,  tho  flow  enters  at  Mach  one  and  accelerates  to  supersonic 
velocities  as  a  portion  of  tho  gas  expands  and  turns  into  the  hole.  The  shock 
at  the  downstreas  lip  of  the  hole  turns  the  expanded  flow  parallel  to  the  solid 
surfaces  and  reduces  the  velocities  to  subsonic  levels.  Tho  pressure  on  the  lip 
-  is  nearly  twice  the  static  pressure  of  the  incooing  stress.  The  flow  acceler¬ 
ates  away  fros  this  region  and  leaves  tho  tube  and  hole  at  supersonic  veloci¬ 
ties.  There  is  a  large  subsonic  region  in  the  upstreas  portion  of  the  hole.  A 
sore  cospleto  description  of  the  flow  and  a  cosparison  with  shock  tube  data  are 
given  in  Reforence  6. 
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Figure  2.  Velocity  vector  and  pressure  contour  plots  for  a  vent  floM 
Nith  a  Mach  nuaber  of  unity  at  the  entrance  plane. 


Figure  3.  Velocity  vector  encl  pressure  contour  plots  for  e  vent  floM 
Nith  •  Meeh  nunber  of  two  et  the  entrence  pisne. 

To  eelctilete  the  blest  field,  the  floN  through  esch  vent  Is  required  at 
each  Instant  of  tlee  during  tube  bloMdown.  Because  the  floM  Is  three- 
dieenslonal,  It  Is  not  practical  to  obtain  the  coeplete  solution  with  a  tran¬ 
sient  calculation.  Fortunately,  the  floM  contains  eany  features  iMhlch  permit  a 
vigorous  simplification  of  the  problem. 

First,  because  of  the  large  volume  of  the  gun  tube,  the  blowdown  process 
takes  on  the  order  of  tens  of  milliseconds,  while  the  three-dimensional  calcula¬ 
tions  Indicate  that  the  flow  In  a  hole  Is  established  In  a  fraction  of  a  milli¬ 
second.  Therefore,  the  latter  can  be  treated  as  quasi-steady  and  only  the  flow 
inside  and  outside  of  the  tube  must  be  considered  as  time-dependent. 
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Steondly,  in  thn  applicatlont  of  intoroot»  tho  flow  it  oithor  oonie  or 
■uportonic  m  it  ontort  tho  broko  ond,  duo  to  tho  vonting,  oxpondo  to  highor 
Koch  nuobors  os  it  trovolt  dOMnstrooo.  Also,  bocouso  of  tho  nigh  tubo 
prossuros,  tho  gos  oxits  ooch  holo  ot  noor  sonic  or  suporsonie  voloeitios  ovor 
oost  of  tho  OMit  plono  oroo  (soo  figuros  2  ond  3).  Cxporioneo  hot  shoMn  thot 
tho  floo  is  rothor  insonsitivo  to  tho  outflou  boundory  condition  ovor  tho 
roMining  subsonic  portion.  Thus,  tho  flou  ot  o  porticulor  holo  locution  is  not 
influoncod  by  ttvonts  occurring  forthor  doMnstroos  or  outside  of  tho  tubo.  It 
doponds  solely  on  tho  conditions  in  tho  tubo  upstroso  of  the  hole. 

Tho  Euler  oqustions  ooy  bo  M'itton  in  consorvstivo  fore  ss 


>9  ♦  55  ♦  5$  ♦  §H  .  0 

at  ax  ay  a*  ” 


Nhoro 


s»/p^P 


F  ■  on/p 


n*/p^P 


H  «  In/p 


IVp+P 


.(E+P)s/p. 


.(E4P)n/p. 


.(E^P)l/p_ 


p  is  tho  density:  o  >  pu,  n  >  pv,  ond  I  ■  pM  sro  tho  oonontuo  conpononts  in  tho 
X,  y,  ond  z  directions,,  rospoctivoly;  u,  v,  end  n  are  the  corresponding  velocity 


components.  P  is  the  pressure  and  E  is  the  total  energy  per  unit  volute  defined 


E  -  pe  ♦  (m«+n«+l«)/2p  (2) 

Mhore  e  is  tho  specific  internsl  energy.  Both  the  propellent  gas  and  tho  air 
are  taken  to  bo  perfect  gasos,  so  the  pressure  is  related  to  the  state  variables 
p  and  e  by  the  expression 

P  -  (y-l)pe  (3) 

Mhore  7  is  the  specific  heat  ratio. 
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Consider  the  result  of  nondimensionaMrlng  th*  Euler  equations  in  the 


folloMing  way: 

P'  ■  P/Pg  »  ■'  ■  »/f^P2P2  »  f''  ■  n//p2P2  »  *'  ■  l/f^P2P2 
P'  ■  P/P2  »  E'  ■  E/P2  .  e*  ■  ep2/P2  >  C  ■  C/P2/P2 
X'  ■  x/D  ,  y'  ■  y/D  ,  2'  ■  z/D  ,  t'  ■  t/p^pJ/D 
where  P2  and  P2  are  the  density  and  pressure  of  the  uniform  flow  at  the  upstream 
plane  of  the  tube  and  D  is  the  vent  diameter.  The  form  of  the  Euler  and  state 
equations  remains  unchanged,  while  the  inflow  boundary  conditions  become 


P'  -  1  (4) 

«•  -  /yM2  (5) 

n*  .  I*  (6) 

E'  -  l/(y-l)  +  yM5/2  (7) 


Since  the  flow  depends  only  on  the  inflow  boundary  conditions,  which  are 
completely  described  by  the  upstream  Hach  number,  M<>,  the  specific  heat  ratio, 
y.  and  the  hole  geometry,  one  solution  with  these  parameters  specified  is  valid 
for  all  upstream  pressures  and  densities.  This  observation  is  central  to  the 
success  of  the  analysis  because,  while  a  wide  range  of  physical  states  are 
encountered  during  blowdown,  only  a  few  three-dimensional  solutions  are  required 
to  describe  them. 

Data  from  the  three-dimensional  solution  are  used  to  obtain  average  values 
of  the  density,  p^.  pressure,  Ph,  the  mass  flux,  phvh*  and  the  axial  component 
of  the  radial  momentum  flux,  p^vhuh,  in  the  exit  plane  of  the  vent  using  the 
following  expressions: 


Ph  ■  (1/Ah)  P'dA  (8) 

Ph  -  (I/Ah)  /a„  P'^IA  (9) 

PHVH  •  (1/Ah)  p'v'dA  (10) 

PHVHUH  ■  (1/AH)  /.  p'v'U'dA  (11) 
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Mh«rt  u  It  tht  vtlocity  cotpontnt  ptralltl  to  tht  tub*  axis  and  v  Is  the  con- 
ponant  paraHal  to  tha  vant  axis.  Tha  Intagritlon  Is  carriad  out  ovar  tha  vant 
axit  araa,  Ah.  Tha  fluxas  ara  calculatad  rathar  than  tha  valocitlas  bacausa 
thay  ara  usad  directly  to  conputa  tha  local  vant  rata  in  tha  intarior  solution 
and  tha  waapon  impulse  (ref  8).  Tha  avaragas  ara  dimansionlass  and  ara  func¬ 
tions  of  tha  paraaatars  that  appaar  in  tha  thraa-dinansional  solution.  Thair 
usa  Mill  ba  dascribad  balow. 
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Figura  4.  Tha  avaragad  functions  conputad  from  tha  thraa-dinansional  solutions. 
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Th«  gas  venting  through  the  brake  Mill  e'ither  be  the  air  ahead  of  the 
projectile— the  precursor  f10N — or  propellant  gas.  Experience  has  shown  that 
the  three-diaensional  solutions  are  not  particularly  sensitive  to  the  value  of 
y,  therefore  the  numerical  results  obtained  in  Reference  8  for  y  «  1.22  were 
used  in  this  study.  The  interior  and  exterior  flows  are  computed  using  the 
appropriate  specific  heat  ratios,  however. 

The  averaged  functions  are  shown  in  Figure  4.  The  hole  geometry  is  charac 
terized  by  the  ratio  of  its  height,  L^,  to  its  diameter,  Dh.  The  height  is 
equal  to  the  tube  wall  thickness.  Note  that  the  momentum  flux  is  negative  for 
both  vent  heights,  especially  the  shorter  one.  This  is  consistent  with  the 
velocity  vector  plots  of  Figures  2  and  3 — the  flow  leaving  the  shorter  vent  is, 
on  balance,  directed  more  upstream. 

THE  INTERIOR  FLOW  FIELD 

The  flow  inside  the  tube  is  calculated  using  the  one'dimensional  Euler 
equations  with  a  source  term  included  to  represent  the  venting  at  the  tube  wall 

a 

ap  am  1  dm 

—  ♦  —  - -  (12 

ax  ax  A  dx 

a 

am  a(m*/p+P)  m  dm 

—  + - - -  (13 

at  ax  pA  dx 

BE  a(m(E+P)/p;  (E-I-P)  dm 

at  ax  pA  dx 

a 

The  vent  term  (l/A)dm/dx  represents  the  mass  of  fluid  per  unit  time  per  unit 
volume  leaving  the  tube  at  x.  A  is  the  bore  area.  The  fluid  is  assumed  to 
leave  at  the  local  velocity  u  in  the  momentum  equation  and  with  the  local 
enthalpy  per  unit  mass  (E^-Pl/p  in  the  energy  equation  (see  Reference  10  for  the 
derivation  of  these  equations). 


Th«  theory  h111  be  coepered  Mith  en  exper-ieentel  brake  having  holes  of  a 
single  diameter  and  a  uniform  spacing,  along  the  tube.  The  vent  area  per 
unit  length  is  sNcOh'/^Sh  where  Nq  is  the  number  of  holes  per  row  (rows  run 
around  the  circumference  of  the  tube).  The  vent  tern  can  then  be  written  in 
dimensional  form  as 

(l/A)dm/dx  «  -PhvH*^(5<C/Sh)(Dh/Db)*  (15) 

where  Dg  is  the  tube  diameter.  P  and  p  are  the  local  values  of  pressure  and 
density  in  the  one-dimensional  solution;  they  appear  since  these  quantities  were 
used  to  nondinensionalize  the  product  phvh  in  the  three-dimensional  solution. 

The  latter  is  also  a  function  of  the  local  Mach  number. 


THE  EXTERIOR  SOLUTION 


The  flow  outside  of  the  tube  is  treated  as  axi symmetric.  The  large  number 
of  vents  typical  of  such  brakes  and  their  symmetrical  placement  around  the  tube 
makes  this  feasible.  The  Euler  equations  take  the  following  form: 
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(16) 


where 
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Since  two  gases  are  present,  the  state  variables  must  be  evaluated  for  a 
aixture.  The  details  are  fully  explained  in  Rererence  11.  The  last  equation  in 
this  set  is  a  species  equatic;i  where  S  represents  the  mass  concentration  of  the 
air  at  a  particular  location. 
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At  th«  nuzzle,  the  one-dlnensionel  equet^ons  ere  solved  sinultaneously  with 


the  axisyenetrie  equations  downstreen.  The  renainder  of  the  tube  is  a  solid 
boundary  except  for  the  vented  region.  Since  the  area  of  each  vent  represents 
only  a  portion  of  the  local  tube  area,  the  averaged  variables  at  the  vent  exit 
have  to  be  adjusted  to  provide  an  appropriate  boundary  condition  for  the  axisym- 
metric  equations. 

The  sketch  in  Figure  5  represents  a  cross  section  through  the  vented  region 
of  the  tube  along  a  plane  parallel  to  the  tube  axis.  Three  rows  of  holes  are 
shown.  The  center  row  has  an  annular  control  volume  drawn  above  it  which 
extends  completely  around  the  circumference  of  the  tube.  The  width  of  the 
control  volume  is  the  hole  spacing,  S^,  and  its  height  is  6.  If  6  is  allowed  to 
approach  zero,  then  the  only  relevant  flux  vector  in  Eq.  (16)  is  6(Q).  The  flow 
adjusts  instantaneously  from  the  smaller  vent  exit  area,  A^,  to  the  larger  tube 
surface  area  A3  >  irlOg  2Lh)Sh/Nc.  Writing  this  balance  out  gives 


PhvhAh  ■  Ps^s^s  (17) 

Pa(As”AH)  *  Ph)Ah  ■  (PS'^S* 

PHVhUhAh  ■  Ps^S^sAs  (19) 

+  Ph)''hAh  ■  (^S  ■*’  PsJ'^SAs  (20) 


The  species  equation  is  not  written  out  since  the  vented  gas  will  either  be  air 
with  83  a  P3  or  propellant  gas  with  S3  >  0. 


Db/2. 


Figure  5.  Control  volume  used  to  construct  the  boundary  conditions 
along  the  vented  portion  of  the  tube. 
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The  first  tsns  on  tho  loft-hsnd  sido  of  tho  sxial  MOMontum  equation  repre- 
sants  soae  average  pressure,  P^,  acting  on  the  area  difference  (As  -  Ah).  Since 
this  pressure  is  not  knoMri,  it  Mill  be  assusied  that  the  gas  undergoes  an 
isentropic  expansion  so  that  this  equation  can  be  replaced  with  the  expression 

Ps/Ph  •  (Ps/Ph)^  (21) 

where  y  is  appropriate  to  the  gas  being  vented.  With  a  little  algebra,  the 
following  results  are  obtained: 


Ms 


1  Lull 

Ah  2  ♦  (y-i)MH*  2  (y-1) 


(22) 


PS 


PH  ts^] 


(23) 


Vs  ■  Vh(PhAh/PsAs)»  Us  ■  Uh  (24) 

where  M  »  u//yP/p.  After  deternining  Ms  by  iteration,  the  conserved  variables 
in  the  Euler  equations  can  be  formed.  The  quantities  at  the  vent  exit  are 
related  to  the  interior  flow  through  the  averaged  functions  given  above. 

Because  the  vent  exit  flow  is  supersonic,  the  exterior  boundary  condition  is 
completely  determined  by  the  local  conditions  in  the  tube. 

Marten's  Total  Variation  Diminishing  scheme  (ref  12)  was  used  in  conjunc¬ 
tion  with  a  time-splitting  algorithm  to  solve  the  Euler  equations.  A  more 
thorough  description  pf  the  numerical  methods  is  given  in  References  6  and  11. 


THE  STARTING  SOLUTION 

The  first  two  configurations  of  Reference  1  will  be  analyzed:  the  bare 
muzzle  case  and  brake  #1.  This  brake  has  16  rows  of  holes  arranged  in  a 
staggered  pattern  with  every  other  row  of  12  holes  rotated  15  degrees  with 
respect  to  the  adjacent  rows.  The  holes  are  drilled  perpendicular  to  the  tube 
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•xls.  Th«  brtk*  adds  11.84  cni  to  tho  longth  of  tho  cannon.  A  aora  complata 
daacription  of  tha  axparlaantal  satup,  including  photographs  of  tha  brakas,  is 
givan  in  Rafaranca  1.  Tha  data  in  Tabla  1  ara  a  iso  takan  fron  that  source. 


TABLE  I.  INITIAL  DATA  FOR  20-MM  CANNON 


Pm  •  287.0  atm 

U  ■  41.7  cm* 

Yp  •  1.25 

Vm  «  1045.0  m/sec 

C  «  0.0389  kg 

Mp  -  22.8 

T|^  ■  143.0  cm 

W  •  0.098  kg 

Dh  ■  0.0579  cm 

Lp  ■  7.5  cm 

Ta  -  1.40 

■  0.0866  cm 

Dg  ■  2.00  cm 

Mm  •  29.0 

■  0.0709  cm 

Pm  and  Vm  ara  tha  naasurad  projectile  base  pressure  and  velocity,  respectively, 
at  the  auz2le  exit  plane.  T|_  is  the  projectile  travel  length,  Lp  is  the  projec¬ 
tile  length,  Og  is  tha  bore  diameter,  U  is  the  chamber  volume,  C  is  the  charge 
mass,  and  W  is  the  projectile  mass,  r  is  the  specific  heat  ratio  and  H  is  the 
molecular  weight  of  tha  gases;  the  subscripts  "a”  and  "p"  refer  to  tha  air  and 
propellant  gas,  respectively. 

The  starting  configurations  are  shown  schematically  in  Figure  6.  The 
overall  length  of  156.3  cm  in  tha  bare  muzzle  case  includes  a  uniform  extension 
of  the  tube  to  account  for  the  chamber  volume.  The  projectile  velocity  is  kept 
constant  throughout  the  calculation. 

Initially,  the  projectile  nose  is  21.8  cm  from  the  nuzzle  with  the  precur¬ 
sor  shock  at  the  muzzle  exit  plane.  This  distance  places  the  precursor  shock  at 
the  correct  position  in  the  first  shadowgraph  to  be  presented  below.  The  state 
of  the  air  ahead  of  the  projectile  corresponds  to  a  shock  Mach  number  based  on 
the  ratio  of  the  projectile  velocity  to  the  acoustic  speed  in  the  stagnant 
environment.  This  method  of  introducing  the  precursor  flow  into  the  calculation 
underestimates  the  mass  of  air  expelled  from  the  weapon  which  is  (T|,-Lp)A  p^. 


it  th«  dtntity  of  tht  tnvironMnt.  In  tho  computation,  tha  21.S**ca  Jir 
column  hat  a  dtntity  of  4.52  p«.  Taking  tha  ratio  of  tha  tMO  mattat  thOMt  that 
38  parcant  aora  air  it  axpallad  by  tha  maapon.  Tha  daficiancy  could  ba 
corraetad  by  uting  a  projactila  trnval-tiaa  curva  from  an  intarnal  ballittict 
coda  to  ganarata  tha  pracurtor  floN  from  thot  start,  but  tha  computation  time 
mould  alto  Incraasa  tignif icantly.  Tha  affact  of  tha  precursor  f1oM  on  tha 
blast  fiald  mill  ba  axplorad  balom.  Idantical  starting  conditions  ware  u«ed 
whan  tha  brake  was  added,  as  indicated  in  Figure  6 


BREECH  PROJECTILE  SHOCK 


BARE  HUZZLE  CONFIGURATION 


IE 


h 


11.84 — 
8.94-1-0.4 


II 


TnOTlIll 


BRAKE 


configuration  uith  brake 


Figure  6.  The  starting  configurations  with  and  without  the  brake. 
All  dimensions  are  in  centimeters. 


Tha  calculation  proceeds  until  the  projectile  base  reaches  the  position  of 
tha  bare  nuzzle  exit  plane.  At  this  instant,  tha  state  of  tha  propellant  gas 
behind  the  projectile  is  computed  from  the  information  in  Table  I.  The  distri¬ 
butions  of  pressure,  density,  and  velocity  are  based  on  tha  Pidduck-Kent 
limiting  solution  (ref  13)  specialized  to  a  perfect  gas.  These  are  given  by 


Q  .  12:11  9  (12:11  9  ♦  (l:Sy!)i/(y-i)du)-i  cas 

“  2y  W  '  Ir  w  'o  '  1-0  '  ' 

X,  -  {\  ♦  U/A)  (26 

X  •  x/X,  (27 

p  .  p.  (l.:,»!3y/(r-i)  (28 

«  .  9  ^2x _ Q.,  w  f.Q 

^  U  (Y-l)  (i-qJ  c  ^  -  q  ^  '‘® 

u  «  XVn  (30 

Q  is  obtMin«d  by  solving  Eq.  (25)  by  itoration.  Xg  is  ths  projsctila  bssc  posi 
tion  which  includas  the  unifore  cheeber  voluee  extension.  It  reeeins  unchanged 
when  the  brake  is  added. 


RESULTS 

A  set  of  five  shadowgraphs  was  obtained  for  each  of  the  brakes  tested  by 
Oi'ilcn  (ref  1)>  but  they  were  not  included  in  that  report.  Soee  of  these  shad¬ 
owgraphs  are  presented  in  this  section — the  rest  may  be  found  in  the  Appendix. 

The  density  contour  plot  in  Figure  7  shows  the  precursor  flow  just  after 
the  projectile  clears  the  muzzle.  As  noted  earlier,  the  starting  configuration 
was  chosen  so  that  the  intersection  the  precursor  shock  with  the  tube  axis 
would  closely  match  that  in  this  shadowgraph.  In  all  of  the  contour  plots  pre¬ 
sented,  the  data  were  dumped  at  the  instant  the  projectile  reached  its  position 
in  the  corresponding  shadowgraph.  Thus,  what  is  being  compared  between  a  con¬ 
tour  plot-shadowgraph  pair  is  the  structure  of  the  remainder  of  the  flow  field. 
The  computed  shock  shape  and  plume  size  compare  reasonably  well  with  the  experi 
ment  in  Figure  7. 


5 


Th*  shadoMgraph  In  Figur*  6  thoMS  that  tha  bUst  M«va  doat  not  hava  tha 
taooth  appaaranea  of  tha  pracuraor  ahock.  Tha  prasaura  contour  plot  In  tha 
figura  also  has  this  faatura.  It  is  dua  to  tha  nonunifora  anvironaant  ganaratad 
bv  tha  pracursor  floM.  This  can  ba  daaonstratad  by  starting  tha  calculation 
with  tha  projactila  nosa  Just  upstraaa  of  tha  auiila,  tharaby  aliainating  tha 
pracursor  floM.  Tha  saooth  blast  aava  and  siaplar  pluaa  structura  of  Figura  9 
rasult.  Including  tha  pracursor  floM  incraasas  tha  coaputation  tiaa  by  60  par- 
cant,  but  it  appaars  to  ba  nacassary. 

Tha  flow  fiald  davalopaant  with  tha  braka  in  placa  is  dapictad  in  Figures 
10  through  12.  In  tha  density  plot  of  Figura  10,  tha  expansion  of  tha  pracursor 
flOM  in  tha  tuba  dua  to  tha  venting  through  tha  braka  is  evident.  It  does  not 
appear  in  Figures  11  and  12  because  tha  contour  levels  Mara  chosen  to  aaphasiza 
tha  structura  of  tha  exterior  floM. 

Except  for  tha  region  where  tha  blast  wave  intersects  tha  tuba  wall,  tha 
gaoaatrias  of  tha  pracursor  and  aain  shocks  show  good  agraaaant  with  those  in 
tha  shadowgraphs.  Tha  weak  shocks  resulting  fro«  tha  interaction  of  tha  brake- 
driven  and  auzzla-drivan  portions  of  tha  flow  are  also  captured  nicely. 
Unfortunately,  tha  danse  propellant  gas  conceals  tha  rich  structura  of  tha 
pluses  in  tha  shadowgraphs. 

Perhaps  the  nost  intriguing  feature  in  the  figures  is  that  the  primary 
thrust  of  the  brake  flow  is  directed  considerably  rearward  even  though  the  holes 
are  drilled  perpendicular  to  the  tube  axis.  This  is  tha  result  of  the  expansion 
of  the  propellant  gas  remaining  in  the  tube  as  it  passes  the  vented  region.  The 
pressure  drop  along  the  tuba  axis,  shown  in  Figure  13,  is  considerable.  Also 
shown  in  this  figure  is  the  pressure  distribution  along  the  outside  surface  of 
the  tube,  that  is,  the  boundary  condition  Ps  in  Eq.  (21).  Ps  depends  directly 
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Figur*  10.  Density  contour  plot  and  shadowgraph  for  the  brake  case, 
t  >  0.213  Msec. 
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Figure  11.  Pressure  contour  plot  end  shadowgreph  for  the  brake  case, 
t  «  0.423  ;«sec. 


20 


I 


on  the  tube  pressure,  so  Its  variation  with  distance  Is  geometrically  similar  to 
that  In  the  tube.  It  Is  well  above  atmospheric  pressure  at  the  upstream  end  of 
the  brake,  but  only  slightly  above  It  at  the  downstream  end.  Relative  to  the 
pressure  level  between  the  plume  boundary  and  the  main  blast  wave,  the  gas  at 
the  upstream  end  of  the  brake  Is  underexpanded,  while  that  at  the  downstream  is 
overexpanded. 


Figure  13.  Pressure  distribution  on  the  tube  axis  through  the  vented  region 
(left)  and  outside  of  the  brake  (right),  t  »  0.546  msec. 

The  velocity  plot  in  Figure  14  shows  the  motion  of  the  propellant  gas 
within  the  plume.  Every  fourth  vector  is  plotted  to  render  some  clarity  to  the 
figure.  The  plume  boundaries  are  indicated  by  the  solid  lines  which  are 
actually  a  plot  of  the  mass  fraction  of  air  with  contour  levels  spaced  closely 
around  0.5.  They  also  appear  in  the  pressure  plots. 


Figure  14.  Velocity  vector  plot  for  the  brake  case,  t  =  0.546  msec. 

The  vectors  indicate  that  the  gas  leaves  the  brake  with  a  slight  upstream 
bias  over  most  of  its  length.  At  the  upstream  end,  the  underexpanded  flow  is 
gradually  turned  nearly  perpendicular  to  the  tube  by  a  weak  compression  shock. 

At  the  downstream  end,  the  overexpanded  flow  is  turned  immediately  upstream  by  a 
compression  shock.  The  shock  system  which  terminates  the  supersonic  flow  in  the 
plume  also  turns  the  flow  upstream.  It  then  undergoes  an  expansion  near  the  top 
of  the  plume  where  the  flow  is  turned  nearly  parallel  to  the  tube  axis.  This 
expansion  is  terminated  by  another  shock  and  the  flow  then  circulates  around  a 

large  vortex  on  the  left  side  of  the  plume.  Note  the  absence  of  a  vortex  on  the 
right  side. 

To  determine  the  influence  of  the  muzzle  flow  on  the  shape  of  the  brake 
plume,  the  calculation  was  repeated  with  the  muzzle  moved  farther  downstream. 

The  pressure  plot  in  Figure  15  shows  that  the  brake  plume  develops  essentially 


as  before. 


Figure  15.  Pressure  contour  plot  for  the  brake  case  with  muzzle  flow  removed, 
t  «  0.546  msec. 

Removing  the  precursor  flow  by  placing  the  projectile  nose  just  upstream  of 
the  brake  at  the  start  of  the  calculation  produced  the  result  in  Figure  16.  The 
plume  shape  is  modified  only  slightly.  It  is  interesting,  however,  that  now  the 
upstream  portion  of  the  blast  wave  more  closely  resembles  the  shadowgraph.  This 
suggests  that  it  may  be  important  to  include  the  full  development  of  the  precur¬ 
sor  flow  from  shot  start  to  obtain  proper  results  in  this  region.  It  should 
also  be  mentioned  that  failure  to  include  the  45-degree  ramp  on  the  tube 
exterior  just  upstream  of  the  experimental  brake  (see  shadowgraph  in  Figure  10) 
might  contribute  to  the  poor  agreement,  although  the  corresponding  step  in  the 
bare  muzzle  case  (see  Figure  7)  does  not  seem  to  influence  either  the  precursor 
or  the  ma^n  shock- development  there. 


Figure  16.  Pressure  contour  plot  for  the  brake  case  with  muzzle  and 
precursor  flows  removed,  t  =  0.54t  msec. 

CONCLUSIONS 

In  general,  the  predicted  results  compare  favorably  with  the  shadowgraphs. 
The  coupling  of  the  interior  and  exterior  flows  using  a  control  volume  model  to 
make  the  transition  appears  to  be  practical.  Perhaps  more  attention  should  be 
paid  to  the  development  of  the  precursor  flow. 
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APPENDIX 


Th«  add-itional  shadowgraphs  referred  to  at  the  beginning  of  the  Results 
Section  are  given  below.  The  data  for  the  first  contour  plot  for  each  con¬ 
figuration  were  dunped  at  the  Instant  the  Intersection  of  the  precursor  shock 
with  the  tube  axis  Matched  that  In  the  corresponding  shadowgraph.  This  was 
necessary  because  the  location  of  the  projectile  In  the  tube  was  not  known.  In 
the  remaining  plots,  the  projectile  posi  on  was  used  as  the  reference,  as 


described  In  the  text. 
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Figure  A-1.  Density  contour  plot  and  shadowgraph  for  the 
bare  muzzle  case,  t  >  0.187  msec. 
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Pressure  contour  plot  and  shadowgraph  for  the 
bare  muzzle  case,  t  «  0.394  msec. 


Pressure  contour  plot  for  the  bare  nuzzle  case,  t  =  0.588  nsec 
corresponding  to  the  shadongraph  in  Figure  A-4. 


Figure  A-4.  ShadoMgraph  for  the  bare  mizzle  case  corresponding 
to  the  pressure  contour  plot  in  Figure  A-3. 


Figure  A-5.  Density  contour  plot  and  shadowgraph  for  the 
brake  case,  t  «  0.315  nsec. 

i 


Figure  A-6.  Pressure  contour  plot  for  the  broke  case,  t  •  0.584  esec 
corresponding  to  the  shedoNgraph  in  Figure  A-7. 


Figure  A-7.  ShadoMgraph  for  the  brake  case  corresponding  to  the 
pressure  contour  plot  in  Figure  A-6. 
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DIRECTOR,  PROCUREMENT  DIRECTORATE  1 

ATTN:  SMCWV-PP 

DIRECTOR,  PRODUCT  ASSURANCE  DIRECTORATE  I 

ATTN:  SMCWV-QA 


NOTE:  PLEASE  NOTIFY  DIRECTOR,  BENET  LABORATORIES,  ATTN:  SMCAT?-CCB-TL ,  OF 
ANY  ADDRESS  CHANGES. 


TECHNICAL  REPORT  EXTERNAL  DISTRIBUTION  LIST 


NO.  OF 
COPIES 

ASST  SEC  OF  THE  ARMY 

RESEARCH  AND  DEVELOPMENT 

ATTN:  DEPT  FOR  SCI  AND  TECH  1 

THE  PENTAGON 

WASHINGTON,  O.C.  20310-0103 
ADMINISTRATOR 

DEFENSE  TECHNICAL  INFO  CENTER 
ATTN:  DTIC-FDAC  12 

CAMERON  STATION 
ALEXANDRIA,  VA  22304-6145 

COMMANDER 
US  ARMY  AROEC 
ATTN:  SMCAR-AEE 

SMCAR-AES,  BLOG.  321 
SMCAR-AET-0,  BLDG.  351N 
SNCAR-CC 
SMCAR-CCP-A 
SMCAR-FSA 
SMCAR-FSM-E 
SMCAR-FSS-0,  BLDG.  94 
SMCAR-lMl-I  (STINFO)  BLOG. 
PICATINNY  ARSENAL,  NJ  07806-5000 

DIRECTOR 

US  ARMY  BALLISTIC  RESEARCH  LABORATORY 
ATTN:  SLCBR-DD-T,  BLDG.  305  1 

ABERDEEN  PROVING  GROUND,  MO  21005-5066 

DIRECTOR 

US  ARMY  MATERIEL  SYSTEMS  ANALYSIS  ACTV 
ATTN:  AMXSY-MP  1 

ABERDEEN  PROVING  GROUND,  MO  21005-5071 

.  COMMANDER 
HQ,  AMCCOM 

ATTN:  AMSMC-IMP-L  1 

ROCK  ISLAND,  IL  61299-6000 


NO.  OF 
COPIES 

COMMANDER 

ROCK  ISLAND  ARSENAL 

ATTN:  SMCRI-ENM  1 

ROCK  ISLAND,  IL  61299-5000 

DIRECTOR 

US  ARMY  INDUSTRIAL  BASE  ENGR  ACTV 
ATTN:  AMXIB-P  1 

ROCK  ISLAND,  IL  61299-7260 

COMMANDER 

US  ARMY  TANK-AUTMV  R&D  COMMAND 
ATTK?  AMSTA-DDL  (TECH  LIB)  1 

WARREN,  MI  48397-5000 

COMMANDER 

US  MILITARY  ACADEMY  1 

ATTN:  DEPARTMENT  OF  MECHANICS 
WEST  POINT,  NY  10996-1792 

US  ARMY  MISSILE  COMMAND 
REDSTONE  SCIENTIFIC  INFO  CTR  2 

ArrN:  DOCUMENTS  SECT,  BLDG.  4484 
REDSTONE  ARSENAL,  AL  35898-5241 

COMMANDER 

US  ARMY  FGN  SCIENCE  AND  TECH  CTR 
ATTN:  DRXST-SD  1 

220  7TH  STREET,  N.E. 

CHARLOTTESVILLE,  VA  22901 

COMMANDER 

US  ARMY  LABCOM 

MATERIALS  TECHNOLOGY  LAB 

ATTN:  SLCMT-IML  (TECH  LIB)  2 

WATERTOWN,  MA  02172-0001 


NOTE;  PLEASE  NOTIFY  COMMANDER,  ARMAMENT  RESEARCH,  DEVELOPMENT,  AND  ENGINEERING 
CENTER,  US  ARMY  AMCCOM,  ATTN:  BENET  LABORATORIES,  SMCAR-CCB-TL, 
WATERVLIET,  NY  12189-4050,  OF  ANY  ADDRESS  CHANGES. 


TECHNICAL  REPORT  EXTERNAL  DISTRIBUTION  LIST  (CONT'D) 


NO.  OF  NO.  OF 

COPIES  COPIES 


COMMANDER 

US  ARMY  LABCOM,  ISA 

ATTN;  SLCIS-IM-TL  1 

2800  POWDER  MILL  ROAD 
ADELPHI,  HD  20783-1145 

COMMANDER 

US  ARMY  RESEARCH  OFFICE 
ATTN;  CHIEF,  IPO 
P.O.  BOX  12211 

RESEARCH  TRIANGLE  PARK,  NC  27709-2211 
DIRECTOR 

US  NAVAL  RESEARCH  LAB 
ATTN:  MATERIALS  SCI  &  TECH  DIVISION 
CODE  26-27  (DOC  LIB) 

WASHINGTON,  O.C.  20375 


COMMANDER 

AIR  FORCE  ARMAMENT  LABORATORY 
ATTN;  AFATL/MN  1 

EQLIN  AFB,  FL  32542-5434 

COMMANDER 

AIR  FORCE  ARMAMENT  LABORATORY 
ATTN:  AFATL/MNF 

1  EQLIN  AFB,  FL  32542-5434  1 

METALS  AND  CERAMICS  INFO  CTR 
BATTELLE  COLUMBUS  DIVISION 
505  KING  AVENUE 

COLUMBUS,  OH  43201-2693  1 

1 

1 


NOTE;  PLEASE  NOTIFY  COMMANDER,  ARMAMENT  RESEARCH,  DEVELOPMENT,  AND  ENGINEERING 
CENTER,  US  ARMY  AMCCOM,  ATTN:  BENET  LABORATORIES,  SMCAR-CCB-TL , 
WATERVLIET,  NY  12189-4050,  OF  ANY  ADDRESS  CHANGES. 


